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The pion production by sigma decay and its relation with chiral symmetry 
restoration in a hot and dense matter are investigated in the framework of the 
Nambu-Jona-Lasinio model. The decay rate for the process a — s- 27r to the lowest 
. order in a I /Nc expansion is calculated as a function of temperature T and chemical 

I potential /i. The thermal and nonthermal enhancements of pions generated by the 

decay before and after the freeze-out are present only in the crossover region of the 
. chiral symmetry transition. The strongest nonthermal enhancement is located in 

I the vicinity of the endpoint of the first-order transition. 
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F; ; 1 Introduction 

> . 

•i-H . 

^ ' It is generally believed that there are two QCD phase transitions in hot and dense nuclear 

. matter One of them is related to the deconfinement process in moving from a hadron 

gas to a quark-gluon plasma, and the other one describes the transition from the chiral 
symmetry breaking phase to the phase in which it is restored. Since the deconfinement 
and restored phases are only intermediate states in relativistic heavy ion collisions which 
are considered to be the way to generate the two transitions in laboratories, the key 
problem is how to extract possible signatures of the phase transitions hidden in the final 
state distributions. 

One way to study the signatures of the chiral transition is to compare the properties 
of hadronic systems with and without undergoing a chirally symmetric phase. In recent 
years many such signatures have been proposed, such as the excess of pions due to a rapid 
thermal and chemical equilibration in the symmetric phase p|, the excess of low energy 
photon pairs by pion annihilation in hot and dense medium 0], the dilepton enhancement 
from leptonic decay of sigma[Q, and the enhancement of the continuum threshold in 
the scalar channel^. Most of the signatures are associated with the changes of sigma 
properties at finite temperatures and densities. In the vacuum there is no obvious reason 
for the presence of the scalar meson sigma in the study of chiral properties. However, there 
is a definite need0 for sigma at finite temperatures and densities in the investigation of 



1 



chiral symmetry restoration. With increasing T and n sigma changes its character from a 
resonance with large mass to a bound state with small mass. This medium dependence of 
sigma mass makes important sense not only in the above mentioned particle yields, but 
also in thermodynamics 0. The contribution from sigma to any thermodynamic function 
can for most case be neglected, because it is heavy, while it plays an important role in 
the region around the critical point of the chiral phase transition, since there m„ = 0. 

From the recent QCD model calculations the chiral phase transition in the chiral 
limit with massless pions is of second order at high temperatures, and of first order at 
high densities. Therefore, there is a tricritical point P in the T — fi plane which separates 
the first- and second-order phase transition. In the real world with nonzero pion mass, 
the second order phase transition becomes a smooth crossover and the tricritical point 
P becomes an endpoint E of a first-order phase transition line. In relativistic heavy ion 
collisions, for the choice of control parameters like colliding energy and impact parameter 
such that the freeze-out of the system occurs near the point E, sigma is one of the most 
numerous species at the freeze-out, since it is nearly massless. These light sigmas at 
the freeze-out will lead to large event-by-event fluctuations of final state pions due to 
the (T-exchange in the pion thermodynamic potential. Another direct signature of this 
critical endpoint E is a. nonthermal excess of pions from the hadronic decay of sigma. 
With the expansion of the hadronic system produced in relativistic heavy ion collisions, 
the in-medium sigma mass rises towards its vacuum value and eventually exceeds the vrvr 
threshold. As the cttttt coupling is large, the decay proceeds rapidly. Since this process 
occurs after freeze-out, the pions generated by it do not have a chance to thermalize. Thus, 
it is expected IP] that the resulting pion spectrum will have a nonthermal enhancement at 
low transverse momentum. 

The present paper is devoted to a systematic and detailed investigation of the above 
idea[3 on excess pions produced by sigma decay in a chiral quark model in the whole 
T — fi plane, without the assumption that the freeze-out occurs at the chiral transition. 
In the chiral limit, we know very well that the sigma decay into two pions starts at the 
phase transition point, and the freeze-out happens later. Therefore, the sigma decay 
before and after the freeze-out will result in thermal and nonthermal enhancements of 
pions, respectively. Our motivation here is to see that in which temperature and density 
region of freeze-out there is a strong nonthermal enhancement. 

We proceed as follows. In Section 2 we calculate the sigma decay rate as a function of 
T and /i to the lowest order in a 1/Nc expansion in the framework of the Nambu-Jona- 
Lasinio model (NJL) at quark level[0], and discuss the importance of the Bose- Einstein 
statistics at high T and the dominance of the first-order transition at high fi. The in- 
medium pion mass drops down with decreasing T and fi. When it approaches to its 
vacuum value, the freeze-out happens. With this definition of freeze-out and the help 
of scaling hydrodynamics to describe the expansion of the hadronic system, we study in 
Section 3 the thermal and nonthermal pion enhancements due to a decay. Finally we give 
our summary and point out open problems that remain. 



2 



2 sigma decay rate at finite temperatures and den- 
sities 



One of the models that enables us to see directly how the dynamical mechanisms of 
chiral symmetry breaking and restoration operate is the NJL model applied to quarks [|10 



Within this model one can obtain the hadronic mass spectrum and the static properties 
of mesons remarkably well. In particular, one can recover the Goldstone mode, and some 
important low-energy properties of current algebra such as the Goldberger-Treiman and 
GellMann-Oakes-Renner relations. Because of the contact interactions between quarks 
that is introduced in this model, there is of course no confinement. A further consequence 
of this feature is that the model is non-renormalizable, and it is necessary to introduce 
a regulator A that serves as a length scale in the problem. The sigma decay rate at 
zero density in this model was firstly studied by Hatsuda and Kunihiro||TT|. In order to 



investigate the thermal and nonthermal pion yields by sigma decay, we recalculate here 
the decay rate in the whole T — ^ plane, and analyze its relation with the first-order 
transition in high density region. 

The two-flavor version of the NJL model is defined through the lagrangian density. 



^NJL 



where only the scalar and pseudoscalar interactions corresponding to a and vr mesons, 
respectively, are considered, ip and ip are the quark fields, the operate r is the SU(2) 
isospin generator, G is the coupling constant with dimension GeV~^, and itlq is the current 
quark mass. 

In an expansion in the inverse number of colors, 1/Ac, the zeroth order (Hartree) 
approximation for quarks together with the first-order (RPA) approximation for mesons 
gives a self-consistent treatment of the quark-meson plasma in the NJL model. In the 
chiral limit the tricritical point P is located at Tp = 79 MeV,^p = 280 MeV in the 
T — II plane. In the case of explicit chiral symmetry breaking the endpoint of the first- 
order transition line is shifted with respect to the tricritical point P towards large fi, its 
position is Te = 23 MeV^^iE = 325 MeV. To the lowest order, namely (l/Nc)^^'^, the 
Feynman diagram for the decay process a ^ 27r is sketched in Fig.l. The partial decay 
rate in the rest frame of sigma is given in terms of the Lorentz-invariant matrix element 
M by 

dT„^2n 1 IpI 



dn 327r2 ml 



\M\' , (2) 



where |p| = y — 'm'i is the pion momentum, mo- and m^r are respectively a- and vr-mass. 
Using an obvious notation for the M-matrix element, one has 

M = g^glA^nn (3) 

with 

A^TTTT = TY^ e*"'"'' f -^-^TrS' (q, iuJn)T^S{q + p, iUn + ^)T^S (q, iuj^ + m^) ■ (4) 
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Here Un = (2n + l)7rT, n = 0, ±1, ±2, ... are fermionic Matsubara frequencies, Tj, = i'-f^r 
is the pseudoscalar interaction vertex. The quark propagator is denoted by S'(q, iun) = 
(7o(zu;„ + yu) - 7 ■ q + rriq) / {{iujn + fi)'^ - E^j with the quark energy Eg = ^Jm'^ + q^. In 
Eq.(^), Tr refers to the trace over color, flavor and spinor indices. The pion-quark-quark 
and sigma- quark- quark coupling strengths in Eq.(^ are determined in the model via|jTO[ 



~2frj. . _ <9n^(A;o,0;T,/i) 
9n {J-^N- ■Q^ \kl=mi 

_2,^ 9n,(A;o,0; T,/i), 



where 11.,^ and 11^ are the standard mesonic polarization functions [|l^ for pion and sigma. 



The dynamically generated quark mass mg(T, /i) and the meson masses rriT^iT,^) and 



mcr{T, fi) are calculated using the usual gap equations [|lOl in the Hartree approximation 
for quarks and RPA approximation for mesons. After evaluation of the Matsubara sum 
in Eq.(l), one hasfTl 



8(q ■ p)2 - {2ml + ■p + mi/2- 2mlE^ 



(ml 



4^2) ((m2 - 2q ■ p)2 - mlE^^ 



(6) 



where frix) = (^1 + e^^"^^ is the Fermi-Dirac distribution function. 

Sigma can decay into neutral and charged pions. By considering the exchange contri- 
bution for a — > 2ttq due to the identical particle effect, the total decay rate can be written 
as 

^a^2n {T, fi) = r^^2no {T, /i) + I^a-^n+TT. {T, /i) 

- ^ ^^^~^^-9l9t\A.UT,l^r ( 1 + 2M^)] , (7) 



71 mi \ 2 



where we have taken into account the Bose-Einstein statistics in terms of the distribution 
function fsix) = {e^/'^ — 1^ for the final state pions |1TT| 



From the comparison of aim coupling strength g^^^T^ defined through the lagrangian 
Li ~ gaTTnO'TTTT and the M-matrix element (|^) in the NJL model, we have 

g^^^{T, /i) = 2M(T, /i) = 2g,glA{T, /i) . (8) 

In the numerical calculations we used the dynamical quark mass = 0.32 GeV, 
the pion decay constant /^r = 0.093 GeV and the pion mass = 0.134 GeV in the 
vacuum (T = /i = 0) as input to determine the 3 parameters in the NJL model, namely 
the coupling constant G = 4.93 GeV~^, the momentum cutoff A = 0.653 GeV and the 
current quark mass m^ = 0.005 GeV. 

We first discuss the temperature dependence of go-Kn and ro-^27r in the case of zero 
baryon density which corresponds to the central region of relativistic heavy ion collisions. 
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The coupling strength (^o-TrTr is about 2 GeV in the vacuum and almost a constant in the 
temperature region T < 0.1 GeV, and then decreases smoothly with increasing tempera- 
ture, as shown in Fig.2. The energy condition for the decay process a ^ 2n to happen 
is rrifj > 2m.^ in the rest frame of a. At the threshold temperature T ~ 0.19 GeV de- 
termined by rriaiT) = 2m^(T), g^^^^T, has a sudden jump from about 0.3/ fm to zero. In 
a wide temperature region, the decay rate (|^) goes up gradually from its vacuum value 
~ 0.45//m to the maximum ~ 0.5/ fm at T ~ 0.17 GeV, then it drops down rapidly but 
continuously and becomes zero at the threshold temperature. The fact that the maximum 
decay rate is not in the vacuum where the mass difference between a and vr has the largest 
value, but at T ~ 0.17 GeV is beyond our expectation. It originates from the important 
Bose-Einstein statistics of the final state pions at high temperatures. This can be seen 
clearly in Fig.2 where we show the decay rates with and without consideration of the 
Bose-Einstein statistics. The rapid decrease of ro-^27r at high temperatures is due to the 
behavior of (^o-tttt and the phase space factor ^ml/i — ml/w?^. 

The density dependence of (/o-tttt and ro-^27r is shown in Fig. 3. Since the threshold 
temperature for a decay decreases with increasing chemical potential, the crucial Bose- 
Einstein statistics effect at high temperatures is gradually washed out as the chemical 
potential goes up. There is almost no difference between with and without such statistics 
when /i > 0.25. Instead of the statistics the first-order chiral transition dominates the 
behavior of (^a-TTTr and T^_,2tt at high densities. When the system passes through the first- 
order transition line from the chirally symmetric phase to the broken phase, all the masses, 
mq,m^ and mg. are discontinuous, and jump up and jumps down. Therefore, 
the mass difference between a and vr has a jump from less than 2mj^ to larger than 2171^^, 
and the a decay starts not at the threshold, but at some value beyond the threshold. 
This means that the decay in the first-order transition region happens suddenly with a 
nonzero rate, the temperature dependence of the decay rate behaves like a step function. 
This is shown in Fig. 3 for = 0.?>?>GeV > fiE- When approaching to the endpoint E 
from the first-order transition side, the jump of the mass difference disappear. The decay 
rate becomes continuous, while the coupling strength has still a jump. 

To see the parameter dependence of the above discussion, we considered another set 
of parameters, G = 5.46 GeV-^,A = 0.632 GeV and mo = 0.0055 GeV. It makes shifts 
of the magnitudes of (^o-tttt and ro.^27r, but does not change their shapes, especially the 
behavior around the threshold temperature or the chiral phase transition point. 

Let's summarize the calculations of the a decay rate. When the chemical poten- 
tial is extremely high, the first-order chiral transition results in an almost temperature- 
independent decay rate. For most of the decay processes in the crossover region, the 
maximum decay rate does not locate in the vacuum, but is very close to the threshold 
point. Only when the chemical potential is around 0.25 GeV which is close to the endpoint 
E, the decay rate increases continuously and not very fast in the beginning, and then it 
behaves like a constant. As we discussed in the introduction, the freeze-out of the pions 
may occur after a decay. According to the decay rate shown in Fig.3, when the decay 
happens in the first-order transition region or most of the places in the crossover region 
in the T — fi plane, a large fraction of cr's may decay into pions before the freeze-out, and 
these pions are thus thermalized before their freeze-out, only those pions generated after 
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the freeze-out contribute to the nonthermal enhancement. However, when the a decay 
happens with n ~ Q.2bGeV , most of the pions may be produced after the freeze-out, and 
therefore, there may be a strong nonthermal enhancement in this case. In the next section 
we discuss the nonthermal enhancement quantitatively. 



3 thermal and nonthermal pion enhancements 

Before we calculate the pion numbers with the a decay rate as a function of temperature 
T and chemical potential /i given in the last section, we need to treat T, /i as time- 
dependent. In the baryon-free (yU = 0) region of relativistic heavy ion collisions, the 
Bjorken's scaling hydrodynamics []T3| is often used to describe the space-time evolution of 
hadronic thermodynamics. The time dependence of temperature in this scenario is given 
by 

T = T,{^jy (9) 

with a = 1/3. Here we have defined td to be the time at which a decay begins and the 
threshold temperature satisfies the mass relation mcr{Td) = 2mT^{Td). We will be inter- 
ested in times t > td- In hydrodynamic models|]T3|, M the initial time corresponding to 



the highest temperature of the fluid is normally taken to be ~ Ifm. Since the sigma decay 
starts in the early stage of the fluid, we choose td = 2/m in our numerical calculations. 

For heavy ion collisions with full nuclear stopping, the lifetime of the baryon-rich 
fluid is much shorter compared with that of the baryon-free fluid, the space evolution 
can be considered as a homogeneous 3-dimensional expansion (|^) with a = 1 instead of 
a = 1/3. The real relativistic heavy ion collisions are between these two limits. We still 
use the expansion mechanism(^ as the time dependence of temperature for any chemical 
potential fi, but the exponent a is now /i-dependent and parameterized by 

a{fi) = I 1 + 2-^) a(0) (10) 

with a;(0) = 1/3. Here the maximum chemical potential \imax corresponding to the full 
stopping limit is chosen as 0.4 GeV in the NJL model. 

With the known decay rate r(T, /i) and the relation between temperature and time 
for any chemical potential, the number of cr's present at time t is determined by 

= -r(t)iv.(t) , (11) 

and is therefore 

iV.(t) = iV,(td)e~^*/^*'^'*' , (12) 
where Na^d) is the maximum number at the beginning time t^, 

= V{td)n^{td) , 
d^p 1 



^'^(t,) = I (2^)3 e..(T,)/T, _ 1 (13) 
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with the a energy e^j = yrril. + and the volume V{td) of the system at time t^- 

Since each sigma yields two pions, the number of pions generated by a decay at time 
t is related to the number of sigma by 

N^{t) = 2{N,{t,)-N,{t)) . (14) 

As the time t ^ oo, N^it oo) = 0, sigma disappears, and the pion number reaches its 
maximum A^7r(t oo) = 2Ncr{td)- 

The conventional definition for freeze-out is that the time between two successful 
scattering events is larger than the collective expansion time scale [p!5[1, 

Tscatt ^ "Texp ■ (1^) 

'Tscatt can be calculated [0 from the thermally averaged vrvr cross section []T^ in the NJL 
model, and Texp = ^/d^uf^ is simply reduced to Texp = t in the scaling hydrodynamics 



Thus the freeze-out time tj and the freeze-out temperature Tj are determined by Tscatt = tf 
and the expansion mechanism (^. With the two temperature scales and Tf we can 
divide the total pions generated by a decay into thermal and nonthermal pions. For 
Td > Tf, the produced pions in the interval Tf < T < T^ have time to thermalize 
before the freeze-out, they lead to an enhancement of thermal pions, while in the interval 
T <Tf the produced pions do not get a chance to thermalize, they result in a nonthermal 
enhancement of pions with low momentum. For T^ < Tf, a decay yields no thermal pions, 
all the resulted pions are nonthermal pions. 

The thermal and nonthermal enhancements are defined as ratios of the numbers of 
the pions produced by a decay before and after the freeze-out to the number of the pions 
produced directly by the freeze-out of the original thermal pion gas, 

^*'^(^) = ' td<t<tf, 

^ntit) = j^d^r^t^) ' tf^^ (16) 

where N^^^{tf) is determined by the freeze-out temperature Tf, the volume of the system 



at time tf, and the pion energy = ym^ + p^, 

Nf^it,) = l-'(«,)nf («/) , 

Obviously the measurable thermal and nonthermal enhancements are rth{tf) and r„t(t — »• 
oo). 



With the in-medium pion momentum p = y ^ — ^l-, the expansion mechanism (^), 
and the pion distribution (p!^), we determine the normalized momentum distribution of 
the nonthermal pions, 

1 dN., e-^l'''^''Tmt-tf)% 
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and the observable mean momentum of the nonthermal pions, 

J dp J^^e J*/^^ r(t)rft 

Since it is difficult to separate out the nonthermal pions experimentally, we calculate 
also the normalized momentum spectra of the total pions including the direct pions, the 
thermal pions from sigma decay, and the nonthermal pions from sigma decay. 



1 dN, 



total ^ 

X 



total 



dp n^-(t^) + 2|^n,(trf) 

t_dnM_f, ^'i:'nt')<it'\\ dnf^{tf,p) , ^t,^ UA^-lL^(''^'^''ruA 



tfnf^{tf)\ )) dp tf "^"^ ^Up^ ^ 



where we have used 



v{tf) - tj ^'^^ 



in the scaling hydrodynamics |T^, and '^"'^ }p^'^^ is defined as 

dnf'-{tf,p) _ p^ 3 



dp 27r2e^-(^/)/^/ - 1 ■ 



(22) 



We now discuss our numerical calculations about the thermal and nonthermal pion 
distributions. Fig. 4 shows the sigma number and pion number scaled by the corresponding 
maximum value, 

rAt)= =l-r^{t), (23) 

NT,[t oo) 

and the thermal and nonthermal enhancements rth{t) and rnt{t) as functions of time 
scaled by the beginning time td in the case of /i = 0.2GeV. The dashed line located at 
the freeze-out time tf separates the left thermal region and the right nonthermal region. 
48% of the (t's decay before the freeze-out, the rest decay into nonthermal pions after the 
freeze-out. The measurable thermal and nonthermal enhancements compared with the 
number of direct pions are 32% and 30%, respectively. 

The measurable thermal, nonthermal and total pion enhancements 

rntit^^)=2^^^rAtf). 

n,{tf)+rUt-^oo) = 2^^^ (24) 
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are shown in Fig. 5 as functions of chemical potential /i. At a first-order chiral transition 
point any particle undergoes a mass jump. At the low end of this jump, rricr < 2mT^, a 
decay is forbidden, while at the high end, mo- > 2m.,^, the decay begins with a nonzero 
rate. This sudden decay was shown in Fig.3. On the other hand, the large a mass at the 
beginning of decay means a strong suppression of the sigma number N„{tii) which dom- 
inates the amplitudes of the enhancements, it drops down exponentially with increasing 
(T-mass. Compared with the decay happen in the crossover region where the decay begins 
with rrifj = 2m.,^, the thermal and nonthermal pion enhancements are both washed away 
by the first-order chiral transition. This is clearly displayed in Fig. 5. The enhancements 
start at the endpoint E of the first-order transition line and develop in the crossover 
region. In a wide chemical potential interval, the total pion enhancement is very strong, 
f'thitf) + rnt{t oo) ~ 0.5. The shapes of the enhancements are mainly controlled by the 
decay rate shown in Fig.3. The nonthermal and total enhancements have their maximum 
values in the vicinity of the endpoint E, this is in agreement with our qualitative analysis 
in the last section. 

The momentum distribution of the nonthermal pions is shown in Fig. 6 for four differ- 
ent chemical potentials. The nonthermal pions are always centralized in a very narrow 
momentum region around p ~ 0.2QGeV. This narrow distribution leads to a sharp peak 
in the momentum spectra for the total pions, shown in Fig. 7. For low chemical potentials 
the nonthermal momenta and the most appropriate thermal momenta almost coincide, 
while for high chemical potentials the nonthermal peak is separated from the thermal 
peak obviously. The thermal peak is always smooth, but the nonthermal peak is always 
sharp. 

The mean momentum of the nonthermal pions is shown in Fig. 8 as a function of 
chemical potential. Unlike the strong chemical potential dependence of the enhancements, 
especially around the endpoint E, the mean momentum is almost a constant, Pnt ~ 
0.26GeV. This value is not so small as the people expected beforeQ]. 

Our study on a decay depends on both the chiral dynamics which determines the decay 
rate, and the expansion mechanism which is reflected in the time evolution of the temper- 
ature. To see which one is dominant in the above calculations, we have chosen different 
expansion parameters. We changed fj.max to control the chemical potential dependence of 
the expansion rate, and the beginning time td of the decay. While the amplitudes of the 
measurable thermal and nonthermal enhancements change remarkably, their shapes, the 
total enhancement and the mean momentum are insensitive to the expansion parameters. 
The chiral symmetry restoration plays a crucial rule in these observable quantities. We 
even let the freeze-out temperature Tf be the decay temperature T^, this means no ther- 
mal pions produced by a decay, the mean momentum in this case decreases by less than 
10% only. 

4 conclusions 

Since the cr-mass goes up and the vr-mass drops down as the hadronic system cools down, 
the more and more large mass difference between a and tt leads to a decay into two pions 
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when the a- mass exceeds the threshold 2m^. These pions resuh in thermal enhancement 
when they are produced before the freeze-out and nonthermal enhancement when they 
are produced after the freeze-out. Based on the two-flavor NJL model we have presented a 
detailed analysis of the pion production at finite temperatures and densities, and discussed 
the relation between the decay and the chiral symmetry restoration, especially the crucial 
effect of the order of the chiral transition on the thermal and nonthermal enhancements. 

We have found that the chiral symmetry breaking and restoration control the pion 
enhancements of the system. At high densities the number of a's at the beginning of their 
decay is strongly suppressed by the large cx-mass, originating from the first-order chiral 
transition. Compared with the large number of direct pions produced by the freeze-out of 
the thermal pion gas, there are almost neither thermal nor nonthermal pion enhancements. 
In practice, the enhancements begin at the endpoint E of the first-order chiral transition 
line and present only in the crossover region with high temperatures and low densities. The 
nonthermal enhancement in this region develops very fast in the early stage of the decay 
and reaches its maximum value in the vicinity of the endpoint E. For a wide chemical 
potential interval, < /i < 0.25 GeV, both the thermal and nonthermal enhancements are 
remarkably large due to the almost constant decay rate during the whole decay process. 
The maximum enhancement is 63% which agrees well with the value 2/3 estimated by 
simple statistics]^. 

While the total enhancement and the structure of the thermal and nonthermal en- 
hancements seem insensitive to the expansion parameters we have chosen in the calcula- 
tions, the amplitudes of the thermal and nonthermal enhancements depend strongly on 
the expansion mechanism. On the other hand, our results are derived in a particular chiral 
model. Although we expect the characteristics of the enhancements to be of more general 
validity, the numerical results are of course model dependent. Therefore, more reliable 
conclusions need more realistic expansion mechanism and more serious QCD model. 
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Figure Captions 



Fig.l: The Feynman diagram for a decay into 27r to the lowest order in 1/Nc expan- 
sion. The sohd hnes denote quarks and antiquarks, and the dashed hnes denote mesons. 

Fig. 2: The decay rate V^^2'k and the couphng strength g^^^^^ as functions of T at = 0. 
The thick and thin sohd hnes indicate the decay rates with and without consideration of 
Bose-Einstein statistics for the final state pions. 

Fig.3: The decay rate Tcr^2-K as a function of T for different chemical potentials. 
The lines marked with a,h,c,d,e correspond to // = 0, 0.1, 0.2, 0.25, 0.3 GeV < /jLe — 
0.325 GeV, respectively, and the line / is with fj, — 0.33 GeV > jiE- 

Fig.4: The sigma number and pion number scaled by the corresponding maximum 
value, and the thermal and nonthermal pion enhancements as functions of time scaled 
by the beginning time of the decay, at ji — 0.2GeV. The dashed line separates the left 
thermal and right nonthermal regions. 

Fig. 5: The measurable thermal, nonthermal and total pion enhancements as functions 

of fjL. 

Fig. 6: The normalized momentum distribution of the nonthermal pions for four dif- 
ferent chemical potentials. 

Fig. 7: The normalized momentum spectra of the total pions for two different chemical 
potentials. 

Fig.8: The mean momentum of the nonthermal pions as a function of 
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